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Abstract 



This report studies a continuous-time version of the well-known Hegselmann-Krause 
model of opinion dynamics with bounded confidence. As the equations of this model 
have discontinuous right-hand side, we study their Krasovskii solutions. We present 
results about existence and completeness of solutions, and asymptotical convergence to 
equilibria featuring a "clusterization" of opinions. The robustness of such equilibria to 
small perturbations is also studied. 

1 Introduction and Preliminaries 

In the study of social dynamics, which has recently attracted much interest from physicists, 
mathematicians and control theorists, modeling the interactions between individuals is the 
core issue. In many successful models of opinion evolution, interactions are assumed to take 
place only if the opinions of the interacting agents are close enough, say, closer than a certain 
threshold: such models are usually called "bounded confidence" models. In this work, we 
study a continuous-time version, introduced in [4], of the "Hegselmann-Krause" model. 

Problem Statement and Contribution 

Let us consider a population of N agents, indexed in a set X = {1, . . . , N}. Each of them has 
a time-dependent real- valued "opinion" Xi(t), which obeys the following dynamics 



Notice that in the above model agent j influences agent i only if \xi — Xj\ < 1, and that 
the function s, which encodes the coupling between the agent opinions and the bounded 
confidence assumption, is discontinuous. Then, the system's right-hand side is discontinuous, 




(1) 



j ex 



where s : R — > K is defined by 
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and this requires to define solutions in suitable sense: in this paper, we use Krasovskii 
solutions, which are defined in the next section. Using Krasovskii solutions to solve a models 
of opinion dynamics with bounded confidence is the main novelty of our work. 

On one hand, the advantage of using Krasovskii solutions is technical, as a complete 
existence and Lyapunov theory is available for such solutions, and this theory allows us to 
obtain interesting results without making the analysis cumbersome. Moreover the set of 
Krasovskii solutions is "large" , so that the results that we state for Krasovskii solutions also 
hold for other types of solutions as Filippov solutions and Caratheodory solutions, in case 
they exist. On the other hand, the convexification of the discontinuities, which is the main 
feature of Krasovskii's definition, can also be seen as a refinement of the original model, 
which smooths out sharp decision thresholds. 

Besides the novelty in the solutions' definition, a distinctive feature of this report is our 
focus on the fundamental properties of the model, in particular (a) average preservation, (b) 
order preservation, (c) contractivity, (d) existence of Lyapunov functions, (e) robustness of 
equilibria with respect to small perturbations. In the literature about networked systems, 
these properties often play an important role in coordination problems: in our work, they are 
instrumental to derive the significant results about solutions, namely existence, completeness, 
and convergence to a clustered configuration. By clustered configuration we mean a state 
configuration such that, for every pair of agents, either their individual states coincide, or 
they differ by more than 1. 

Related Works 

Recently, a vast literature (cf. the surveys [7, 15]) has been produced about the evolution of 
opinions in social dynamics, and the interaction rules which shape such evolution. Scholars 
have remarked that, in spite of significant social forces towards homogenization and con- 
sensus, disagreement persists in societies. A convincing explanation for this phenomenon is 
based on "bounded confidence" : interactions arc limited to be effective only between indi- 
viduals whose opinions are already close enough. A celebrated bounded confidence model is 
due to Hegselmann and Krause [14]. Similarly to most other models, the model consists of 
a discrete-time dynamical system: discrete time allows for immediate computer simulations 
and other analysis advantages, which are exploited for instance in [3], but it entails the draw- 
back of assuming synchrony among the opinions updates. For this reason, it seems worth to 
consider continuous-time bounded confidence models, in which a fixed time schedule is not 
required. This idea was already developed in [4], and indeed our work is strictly related to 
the one presented in [4, Section 2] and in the auxiliary report [2]. In that pair of papers, the 
authors study the same bounded confidence model as here, but they assume a more restric- 
tive (stronger) definition of solutions: this work provides an extension of their results. The 
relationship of our results with the mentioned pair of papers is further discussed in Section 3. 

Finally, we note that a conference version of this work has appeared in the Proceedings 
of the 18th World Congress of the International Federation of Automatic Control as [9] and 
that a related paper focusing on the role of the discontinuity of s is going to appear as [8] . 
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Notation and Preliminaries 

Graphs 

In this paper, we shall make use of some notions from graph theory, and in particular from 
algebraic graph theory. Indeed, graph theory provides an effective tool to model interactions 
between agents and its use is becoming common both in engineering [5, 16] and in economics 
and social sciences [11]. A (weighted) graph G is a triple (V, E, A) where V is a finite set 
of vertices or nodes, E C V x V is a set of edges and the adjacency matrix A is a matrix 
of weights, such that for any u, v G V, A uv > only if (it, v) S E. The Laplacian matrix 
of G is defined as L uv = — A uv when u ^ v and L uu = J2 v gv A U v If G E, then v is 

said to be a neighbor of u in G. A path (of length V) from u to v in G is an ordered list of 
edges (ei, . . . , ej) in the form ((u, wi), (u>i, u>2), (1U2, 103), . . . , (u>i-i, «))■ Two nodes u,v G 
are said to be connected if there exists a path from u to v, and disconnected otherwise. A 
graph is said to be connected if every two nodes are connected, and disconnected otherwise. 
A graph is said to be symmetric when (it, v) G E implies (v, u) G E and the matrix A is 
symmetric. In a symmetric graph, being neighbors is an equivalence relation between nodes: 
the corresponding equivalence classes are said to be the connected components of the graph. 

Solutions to ODEs 

As already remarked, in order to deal with possibly discontinuous ODEs, we need to take 
different notions of solutions into consideration. We provide the definitions of classical, 
Carathcodory and Krasovskii solutions: the interested reader can find more background 
information in classical books as [12] or in the recent tutorial [10]. Let us consider the 
differential equation 



where x : K -> g : R N -> R N , and x G R N . 

A classical solution to (2) on an interval I C M containing to, is a map 4> : I — > M. N such 
that 

1. (j> is differentiable on /, 

2. <j>(t ) = x, 

3. <j>(t) = g{4>{t)) for all t G I. 

A Carathcodory solution to (2) on an interval I C M. containing to j is a map 4> : I — >• R w 




(2) 



such that 



1. (f> is absolutely continuous on /, 



2. 0(t„) = 5, 



3. <p(t) = g(c/>(t)) for almost every t e J. 
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Equivalently, a Caratheodory solution to (2) is a solution to the integral equation 

x(t) = x + / g(x(s))ds. 
J t 

A Krasovskii solution to (2) on an interval / C K containing to, is a map <j> : I —¥ M. N 
such that 

1. <j) is absolutely continuous on /, 

2. 0(to) = x, 

3. </>(i) <E K.g((j)(t)) for almost every is/, where 

£ff(aO = P| co{{g(y) ■ V such that \\x - y\\ < 5}) 

S>0 

and given a set A, by co(A) we denote the closed convex hull of A. 

From the above definitions, it is clear that classical solutions are Caratheodory solutions 
and, in turns, Caratheodory solutions are Krasovskii solutions. Note also that Caratheodory 
solutions coincide with solutions in the classical sense when g is continuous. 

2 Analysis and Results 

This section contains our results about the dynamics (1). After providing a graph-theoretical 
interpretation of the system, we prove several preliminary properties of Krasowskii solutions, 
including completeness and order preservation. Then, we take advantage of these results to 
prove convergence to clustered configurations. Finally, we define the robustness of clusters 
to small perturbations, and provide a necessary and sufficient condition. 

2.1 Interaction Graphs and Basic Properties 

It is useful and suggestive to rewrite system (1) as a dynamic over a suitable state-dependent 
weighted graph, which represents the coupling between the opinions of different agents. In 
such a graph the agents are the nodes, and the opinions of two agents depend on each other 
whenever the agents are neighbors in the graph. By the way system (1) is defined, such 
interaction graph depends on the opinion states via the function s. More precisely, for any 
x € R N we define an interaction graph G(x) — (I, £ (x), A(x)) where the edge set is 

£{x) = i,j el : \xi - Xj\ < 1} , 

that is, S £{x) if and only if s(xj 
by 

, / s \ six 



Xi) > 0, and the adjacency matrix A(x) is defined 
if j =i 
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that is, A(x)ij = 1 if and only if \x% — Xj\ < 1 and j i. The Laplacian matrix L(x) 
associated to Q(x) is then given by 

t, \ [s{xi -Xj) if j ^ i ■ <y 

\Lk±i s ( x k - Xi) it 3=1 

In order to deal with the discontinuity, it is also useful to identify "border" configurations 
by the following definitions of border edge set 

d£(x) = {(i, j), i, j <E I : \xi — xj = 1} . 

and graph Q(x) = (I, £(x),A) , with £{x) = £(x) U d£(x) and A{x)ij = 1 if and only if 
|xi — Xj < 1 and j 7^ z. 

Remark 1 (Symmetry and translation invariancc). We remark that the graphs Q and Q 
are symmetric and invariant with respect to the translation x + al, where a G M and 
1 = (1, 1) T , i.e. = Q{x + al) and §(x) = Q(x + al). 

As we said, the graphs introduced above are interaction graphs in the following sense: if 
two nodes are disconnected, they can not influence each other opinions. 

With the above notation, system (1) can be written as 

x = —L(x)x, (3) 
being L(x) the Laplacian matrix of the state-dependent graph Q{x) and 

(— L(x)x)i = s(xj - Xi)(xj - Xi) 

the components of the right-hand side. As the differential equation (3) has a discontinuous 
right-hand side, we consider Krasovskii solutions to (1), which we characterize as follows. 
For any H C d£(x) we let L H (x) be the Laplacian matrix associated to the graph Q H (x) 
with edges £(x) U H, and correspondingly 

(-L H (x)x) i = J2 Mt)-x t (t)). 

j:(i,j)eS(x)UH 

By the definition, it is clear that a Krasovskii solution to (1) satisfies at almost every time 
the inclusion 

x e co({-L H {x)x : Hcd£(x)}), 

or equivalently the inclusion 

±e{- ^2 &hL h {x)x : a H > 0, ^ ax = lj- 

Hcd£(x) KcdS(x) 

Namely, for a given Krasovskii solution </)(•), 

<j>(t) = - ^2 a^(t)L ff (0(t))0(i) for almost every t, 

Hcd£(<t>(t)) 
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where the time-dependent coefficients a H depend on the solution (/>(•) itself. 

Using this graph-theoretical characterization, we now prove some basic properties of 
Krasovskii solutions to (1). 

Proposition 1 (Basic properties of solutions). Let x(-) be a Krasovskii solution to (1), on 
its domain of definition. 

(i) (Existence). For any initial condition x € WL X , there exists a local Krasovskii solution 
to (1). 

(it) (Order preservation). For any i,j 6 I, if Xi{t\) < Xj{t\), then Xi{t 2 ) < Xj{t 2 ), for any 
t 2 >h. 

(Hi) (Contractivity). For any t 2 > ti, cd({xi(t 2 )} ieI ) C co({xi(ti)}i e x). 
(iv) (Completeness). The solution x(-) is complete. 

(v) (Average preservation). Let x avc (t) = N^ 1 Y^iLi x %if)- Then x avc (t) = x ave (0), for 
t > 0; 

Proof. In the proof, the following notation will be useful. For every i £ X, and every x € M. N , 
we let 

Afi(x) := {k el : \x l -x k \<l} 1 
and for any H C d£(x), we let 

N? = {k G T : (i, k)eH}. 
Clearly, Nf C dAfi(x) := {k E 1 : \xi - x k \ = 1} . With this notation, 

and 

{-L H {x)x) i = (xj(t)- Xi (t)). 
jeA^(x) u Nf 

We are now ready to prove our statements. 

i) Since the right-hand side of (1) is locally essentially bounded, local existence of a 
Krasovskii solution is guaranteed (see for instance [13]). 

ii) To prove the claim, we study the dynamics of the difference between Xj and a;,. By 
continuity of the solutions, we can assume with no loss of generality that Xj and Xi are 
close, for instance that Xj — Xi < 1. For brevity, in the following we omit the explicit 



G 



dependence of A/i on x. For almost every time t, we have 



— (x j -x i ) = aH 

Hed£(x) 



E ( Xh ~ x ^ ~ E ( Xh ~ x ^ 

heMi u nF 



heNj U Nf 



= E aff E U Xh ~ x ^ ~ ( Xh ~ a?i )) + E (( Xh ~ x i) ~ ( Xh ~ Xi ^ 
Hed£{x) [hetfi nAfj fteMniVf 

+ Y i( X h - x j) - ( X h - x i)) - E ( Xfl ~ X ^ 

heNf* n Mj he{Ni u nP)\{M } u n h ) 



+ Y ( X h-Xj) 

he(Nj u N?)\(Ni UN?) 



Hed£(x) 



E -( X j~ X i)+ E {~ X 3+ X i)+ E {~ X j+ X i) 

heMi n Mj heM n n h heNf 1 n M } 



E 



(x h - xCj + 



E 



(X h -Xj) 



he(M U Nf )\(JVj U 2V/ 1 ) 



fcG(A/} U Nf)\(Af, U JVf ) 



= -|Mn^f|(s i -s i )+ E as -(|MnJ\rf | + |JVfnJ\y}|)(s J --x i ) 

HC5£(x) 

E ~ x «) + E ~ x i) 

h&{Ni u N?)\(Nj u JVf ) fee(JVj u Nf)\(Ni u JVf ) 

Since if /i G (A/iUiVf) \ (N?Uj\fj), then x h - X; < 0, whereas if ft G (NjUNf) \ 
(NfllM), then x h - Xj > 0, and since |M n iVf I < |M|, we get that 



j{xj - xi) > nAfj\ + \Mi\ + \AfjDixj - xi). 

The obtained inequality ensures that Xj — Xi can not reach zero in finite time, and yields 
our claim. 

iii) To prove the claim we show that the leftmost agent can only move to its right. To this 
goal, we need a recall the proof of statement (ii) . While our argument shows that strict 
inequalities between agents' states are preserved by the dynamics, we have to remark that 
equalities are not. It is not in general true that if Xi{t{) = Xj(t\), then Xifa) = Xjfo) 
for any t2 > t\. Indeed, we can observe that if Xj(ii) = Xj(t\), then X;(ii) and Xj{t\) 
have to satisfy to the same differential inclusion, but need not to be equal. However, it 
can be proven that it is always possible, given a solution x(-), to sort the states so that 
(t) < Xi 2 (£)<...< Xi N (t), for every t. Note that this mapping i/,) : {1, . . . , N} — > X 
depends on the solution and needs not to be unique. Nevertheless, it allows us to define 
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x m in(t) '■= Xi x {t), and x max (t) := Xi N (t). This fact is useful because it allows us to 
observe that 

Xi(t) - x min {t) > 

for every t and every i £ I and then, for almost every time t, ^x m \ n {t) G [0, +00). 
Repeating an analogous argument for x max implies the claim. 

iv) Claim (iii) ensures that solutions are bounded. By standard arguments, this is enough 
to guarantee that local solutions can be extended for all t > 0. 

v) For every x e M. N , every H C d£(x) and every i, j E 1, it holds that j £ Mi{x) U 

if and only if i € A/} (x) UNj, that is, the graph Q H (x) is symmetric. This key remark 
allows us to argue that for almost every time t, 



=7V_1 E E E (xi(t)-x t (t)) 

iel Hcd£(x(t)) jeA/"i(x)UJVf 

-N- 1 J2 a H(t)J2 E (x j (t)-x t (t))=0. 

Hcd£(x(t)) i£Zj€Mi(x)UNf 



This ensures x avc (t) = x aV c(0) for every t > 0. 



□ 



2.2 Convergence 

We are now ready to prove convergence to a configuration in which agents are separated into 
clusters of agents which share the same opinion. We first recall that a point x is said to be 
a Krasovskii equilibrium of (1) if the function x(t) = x is a Krasovskii solution to (1), i.e. 
G m{{-L H {Z)x : H C d£(x)}). 

Theorem 2 (Convergence of DHK). The set of Krasovskii equilibria of (1) is 

F = ^x € R w : for every (i,j) Elxl, either x% = Xj or \x% — Xj\ > l| 

anrf if x(-) is a Krasovskii solution to (1), then x(t) converges to a point x* G F as t — > +00. 

Proof. The proof is in three steps: we first describe the set of equilibria, then prove conver- 
gence to this set, and finally prove convergence to one equilibrium. 

i) It is clear that every point in F is an equilibrium. To prove that there are no other 
equilibria, we proceed as follows. Without loss of generality we can sort the components 
of x so that Xi 1 < ... < Xi N . For a vector v £ co({— L H (x)Z : H C d£{x)}) to be equal 
to zero, it is necessary that = 0. But since Xk —Xi t > for every k € I, it is necessary 
that Xj — Xi 1 E {0} U[l, +00), for every j € X. Repeating this reasoning for £2, . . ., we 
have that the set of equilibria actually coincides with F. 
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ii) We define the Lyapunov function V(x) = \ ^2i e x x i ana - compute, using the symmetry 
of the graph Q (x) as done in [6] , 

frVW)) =S^(*)i*(t) 

=$><(*)( £ fo(*) -**(*)) + E «h X) 

E fe(*)-^(*)) 2 -^ X a » X (^-w-^(*)) 2 <o. 

(ij)G£(x) Hcd£(x) (i,j)eH 

Since the inequality is strict if x(t) ^ P, and P is closed and weakly invariant, we can 
apply a LaSalle invariance principle [1, Theorem 3] to conclude convergence to the set P. 

iii) We observe that the set F is the union of a finite number of sets Fp, where P = 
{Pi, . . . ,Pfe} is a partition of I in 1 < k < N subsets, and 

Fp = {x G H N : Vi,j gl, if 3ft s.t. i, j e Ph, then xi = Xj,else \xi — Xj\ > 1} . 

As the sets Fp C P are closed and disjoint, each solution converges towards one of 
them. Without loss of generality, we relabel the states so that, for the solution at hand, 
Xi(t) < • • ■ < xjv(i) for every t > 0. When k = 1, the only partition is the trivial one, 
corresponding to equilibria in which the states of all agents coincide. In this case, average 
preservation implies that Xi(t) — > x avc (0) for all i as t — > oo. When fc = 2, there exists a € 
{1, . . . , iV} such that for every x € Fp it holds that Xi — x a for every i < a, Xi = x a+ \ for 
every i > a, and x a +i—x a > 1. Let T a = inf {£ > : x a +i (t) — X a {t) > 1} . If T a < +oo, 
then there is disconnection at finite time and Xi(t) — > ^^2j< a Xj{T) if i < a whereas 
— > jv^a X)j> a x j(T) otherwise. If instead T = oo, then x a +i(t) — x a (t) — > 1 _ as 
t — > +oo. By the average preservation, we argue that Xi(t) — > x avc (0) — if i < a and 
Xi(t) —> a; avc (0) + 4t otherwise. As the argument can be extended to fc > 3 by defining 
k — 1 appropriate disconnection times, we conclude that every solution converges to a 
point in P. □ 

The set of equilibria P in Theorem 2 has the following feature: its points are such that the 
agent opinions either coincide or their distance is larger than 1. Equivalently, the opinions 
of two agents are equal if and only if they are connected in the limit interaction graph. 
Following the opinion dynamics literature, we refer to such groups of agents as clusters, 
and to the corresponding values as cluster values or cluster points. More formally, one can 
consider for a given x G P„ , the map I 9 i i— > Xi G M: the image of such map consists of the 
cluster values, and the clusters are the preimages of the cluster values. The size of a cluster 
is its cardinality. 

Remark 2. (Weak and strong equilibria) According to the definition of Krasovskii equilib- 
rium, Krasovskii solutions which have Krasovskii equilibria as initial conditions may leave 
the equilibria. For example, if N = 2, x = (1,0) G P, there are two Krasovskii solutions 
issuing from x: x 1 ^) = (1,0) and x 2 (t) = (1/2 + l/2e~ 2t , 1/2 - l/2e~ 2t ). In other words, 
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the set F is weakly invariant but not strongly invariant. A subset of F which is strongly 
invariant is F = |a; G R w : for every elxl, either Xi = x j or |xj — Xj\ > l|. 

2.3 Robustness 

In [4], motivated by explaining simulation results about (1), the authors provide a robustness 
analysis for the equilibria of a suitable weighted version of the model. Inspired by this 
approach, we propose a similar definition of robustness for the equilibria of the original 
system (1). Loosely speaking, an equilibrium is said to be robust if the perturbation due to 
adding one agent is not able to make two clusters merge. Equivalently, an equilibrium is said 
to be robust if, after adding one agent, there is no solution whose limit equilibrium point has 
a smaller number of clusters. A more formal definition can be given as follows. 

Definition 1 (Robust Equilibrium). Let x* £ F and xq € K> and consider a Krasovskii 
solution to (1), x(-), such that x(0) = x — (x*, . . . x* N , xq). Let x = lim t _ ) . 00 x(t). 

If, for any xq G M. and any complete Krasovskii solution, the number of clusters of x is 
not smaller than the number of clusters in x*, then x* is said to be robust. 

The following is our main robustness result. 

Theorem 3 (Robustness conditions). Let x* £ F and, when considering any pair of sub- 
sequent clusters in x* , denote their values as xa and xb and their sizes as ua and Hb, 
with ha < Tib. The equilibrium x* is robust if and only if, for every pair of clusters, either 
n-A = n B and \xa — xb\ > 2, or ua < ns and 



1 + — ) (1 + ? ) 

n B J \ n A + n B , 



it 



where the negative number t* n is such that 



e ' n A> n B H 1 H e n A< n B = 1. (4) 

nA + nB nB \ nA + tib J 

The proof of this result is postponed to the Appendix: here we briefly discuss its meaning 
for large populations of agents. 

Corollary 4 (Robustness in large populations). As nA — > oo the necessary and sufficient 
condition in Theorem 3 degenerates into 

\XA ~ X B \ > H 

n B 



for any pair of clusters A and B. 
Proof. From (4) we deduce that 



nA + n B 

and then 



~ ( tn -^ )S n A + n B + l' 
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Consequently, ^1 + ^1 + e tn A>™B is asymptotically equivalent to 1 + ^ as 

?ia — )■ oo. □ 

This corollary states that when clusters are very large, the simple condition \xa — %b \ > 
1 + is necessary and sufficient for robustness against small perturbations. The interest for 
robust equilibria is motivated by the following intuition. Robust equilibria are more suitable 
to be limit points of "real" opinion dynamics system, which would be subject to various 
uncertainties and disturbances. Furthermore, as noted in [4], simulated solutions typically 
converge to robust equilibria. 



3 Krasovskii and Caratheodory Solutions 

This section is devoted to compare the properties of Krasovskii solutions with those of 
Caratheodory solutions, and especially with those of the solutions considered in [4]. 



3.1 Sliding Mode Solutions 

As a consequence of their definitions, the set of Krasovskii solutions may be larger than the 
set of solutions intended in a Caratheodory sense. We now provide an example of a solution 
sliding on a discontinuity surface, proving that there are Krasovskii solutions to (1) which 
are not Caratheodory solutions. 

Example 1 (Sliding mode). Let N = 3 and consider a configuration x in which 1 > X2~ x\ > 
and X3 — X2 = 1. Then, x is on a discontinuity surface due to the disconnection between 
agents 2 and 3. Then, for almost every time 



x € { a 



x 2 - X\ 




X2 


- Xi 


1 + X\ — X2 


+ (l-a) 


X\ 


- X 2 


-1 










a e [0, 1] 



Since the normal vector to the discontinuity plane is vj_ = [0, —1, 1], we have that 



v± ■ x 



-2a + X2 — x\ 



is equal to zero if a = \{x2 — x\). Namely, the Krasovskii solution corresponding to such a 
does not exit the discontinuity plane a; 3 — X2 = 1 at time 0, but it slides on it. The sliding 
solution takes into account the fact that opinions 2:3 and X2 may remain for a while at the 
threshold distance before reaching an equilibrium configuration. □ 

It is an open question whether sliding mode solutions can be attractive for the dynamics. 
However, we know from [4] that a unique complete Caratheodory solution exists for almost 
every initial condition. This implies that the set of initial conditions such that the correspond- 
ing solutions converge to a sliding mode has measure zero, because solutions corresponding 
to those initial conditions would not be complete. 
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3.2 Krasovskii and Proper Solutions 



In [4] and [2], the authors consider a carefully defined subset of Caratheodory solutions to (1): 
they call proper solution any Caratheodory solution x(t) corresponding to an initial condition 
Xq (called proper initial condition) such that 

a) x(t) is the unique Caratheodory solution to (1) with initial condition xq defined on 

[0,+oo); 

b) the subset of [0, +oo) where x(t) is not differentiable is at most countable, and has no 
accumulation points; 

c) if Xi(t) = Xj(t) for some t, then Xi(t') = Xj{t') for all t 1 > t. 

Moreover they prove that almost all x € M. N are proper initial conditions, and that 
proper solutions are contractive in the sense of Proposition 1, preserve the average of states 
and converge to clusters. Our analysis has shown that the most significant properties of 
proper solutions also hold in the larger set of Krasovskii solutions. Nevertheless, there are a 
few significant differences, which we detail in the following list. 

• Existence. For any point xq in WL N , there is a Krasovskii solution x(t) such that 
x(0) = Xo- Instead, there are points, which belong to a certain set P of measure zero, 
such that Caratheodory solutions starting at these points may either not exist or not 
be unique, so that proper solutions may not exist. The set P includes points on the 
discontinuity surfaces, such that x% = Xj for some i ^ j. Note that Krasovskii solutions 
include sliding mode solutions as the one in Example 1: such solutions belong for a 
positive duration of time to a discontinuity surface. 

• Uniqueness. Proper solutions are unique by definition, whereas Krasovskii solutions 
are in general not unique. Note, however, that the results obtained in this note, and in 
particular convergence to a clustered configuration, hold for every Krasovskii solution. 

• Regularity. Krasovskii solutions are differentiable almost everywhere: proper solu- 
tions, by definition, are differentiable out of countable set with no accumulation point. 

• Order preservation. Proper solutions preserve both inequalities and equalities be- 
tween states, while we have remarked that Krasovskii solutions may not preserve equal- 
ities at discontinuities. 

• Connectivity. Along proper solutions, the number of connected components in Q(x(t)) 
is nondecreasing in time. Similarly, along Krasovskii solutions the number of connected 
components in Q(x(t)) is nondecreasing in time. 

• Robustness. We have defined an equilibrium to be robust if the addition of one 
perturbing agent does not make two clusters merge. In [4], robustness is defined for a 
suitable extension of the model, which provides agents with weights: an equilibrium is 
robust if the addition of one perturbing agent of arbitrary small weight does not make 
two clusters merge. Our definition avoids defining this auxiliary weighted system, and 
the small-weight limit is replaced by a limit in the size of the clusters. The resulting 
analysis provides a necessary and sufficient condition, which takes the same simple and 
intuitive form as the condition in [4]. 
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The main drawback of the approach taken in [4] is difficulty in studying existence and 
continuation properties of those solutions. Instead, Krasovskii solutions are easier to deal 
with, as far as existence and continuation properties are considered. Moreover, results about 
proper solutions can be a posteriori obtained as particular cases of the more general results 
on Krasovskii solutions, which also include solutions starting at "problematic points" which 
may not admit proper solutions starting from them. For these reasons, we believe that 
Krasovskii solutions can be a useful tool in opinion dynamics, whenever the model involves 
discontinuities. 
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A Proof of Theorem 3 

Without loss of generality, we restrict our attention to one pair of clusters with values xa, 
xb and sizes ha, tlb, and a perturbing agent with value xq € (xa,xb). We also assume that 
tia < tib and that the agents in each cluster preserve the equality between their states. 1 
Thanks to this assumption, we can limit ourselves to consider the following system of three 
equations, 

' ±b = s(x B - xa)(xq - xb) 

xo = ~n A s(x - x A ) (x - x A ) + n B s(x B - x )(x B - x ) (5) 
± A = s(x A - x )(x - x A ). 

Then, by defining x = xq — xa and y = ig — xq, we are left to study the following bidimen- 
sional (discontinuous) system, 

(x = -(n A + l)s(x)x + n B s{y)y ^ 
\y = n A s(x)x - (n B + l)s(y)y, 

when the initial condition is such that (x(0), y(0)) € (0, 1) x (0, 1). The case of the two original 
clusters A, B merging is equivalent to system (6) converging to the origin. Then, the core of 
our analysis consists in studying for this system the region of attraction of the origin. 

Lemma 5. System (6) has a complete Krasovskii solution converging to the origin if and 
only if the initial condition belongs to the region 1Z delimited by the positive x and y axes, 
the lines {x = 1} and {y = 1} and the branch of the curve 

-{n A +n B + l)t + ^ + 1 j e - t) 1 e -(„ A+ n B +l)t + ^ + 1 y 



n A + n B n A + n B n A + n B UA + n B 

when t* nB < t < and t* nA nB is defined as in (4) 



1 The latter assumption is in general a restriction, as we know from the proof of Theorem 2 that equalities 
between states are not always preserved along Krasovskii solutions at discontinuities. We will show later that 
in this case the assumption entails no loss of generality. 
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Proof. As long as (x(t), y(t)) £ (0, 1) x (0, 1), system (6) reduces to the linear system 



ix = -(n A + l)x + n B y ^ 
\y = n A x - (n B + l)y. 

Being system (7) asymptotically stable, we argue that all solutions, which do not leave 
the unit square, converge to zero. Conversely, outside the unit square the system is not 
asymptotically stable: for instance, if x > 1 and y < 1, we have 

\x = n B y 

\y = -{n B + 

Then, we only need to find which solutions leave the unite square. To this goal, continuity of 
solutions leads us to consider the system in the limit for x or y approaching 1 (from below). 
When y — > 1~ , system (7) becomes 

{x = — {n A + 1)£ + n B 
y = n A x - (n B + 1). 

Since n A x — [n B + 1) < n A — (n B + 1) < 0, solutions can not reach the discontinuity {y = 1}. 
When x — > 1~, system (7) becomes 

\x = -{n A + 1) + n B y 
\y = n A - (n B + l)y, 

implying that solutions may cross the discontinuity {x = 1} if and only if y > . 

Let us then consider the solution to (7), passing at time t = by the point fl, 
which can be written in closed form as 




n A + n B 



-e 



-(n A +n B +l)t _|_ M _)_ 1 ] 

\ n A + n B J 



1 e -( nA +n B +l)t , n A ( 1 1 \- t 



n A + n B n B \ n A + n B 

Note that, by definition, y(t) = 1 if t = t* . As t* A ng is finite, the limited region TZ is 
well defined in the statement of the theorem. By uniqueness of the solutions to (7), we argue 
that a solution to (7) may leave the unit square if and only its initial condition is "above" 
the solution (x,y), i.e., outside the region TZ. □ 

The following corollary derives from Lemma 5 a necessary and sufficient condition for 
solutions to (6) to converge to zero, in terms of the sum a;(0) + y(0). 

Corollary 6. Let d = ( 1 + — ) ( 1 H t — ) e _ *"^'"s . If the initial condition is such that 

J \ ns ) \ nA+riB J J 

x(0) + y(0) > d, then there exists no solution to (6) which converges to zero. Conversely, for 
any d < d, one can find a pair (x(0),y(0)) such that x(0) + y(0) = d and there is a solution 
originating from (x(0),y(0)) which converges to zero. 
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Proof. Let us then consider again the solution to (7) passing by (l, j ■ Note that 

x(t) + y(t)=U + — ) fl + — J— ) 
V n B J \ UA + riB J 

This implies that 

max {x(*) + y(t) : t G [t* n ^ nB , 0]} = x(C A ,„ B ) + v(C»b ) = d ~ 

and 

min{x(t)+y(t) : i e [< Ai „ B , 0]} = £ (0) + y(0) = 1 + ^±1. 

By Lemma 5 and this remark, we deduce the following fact. Being d > fixed, we can find 
(x(0), y(0)) such that x(0) + y(0) = d and (x(0), y(0)) G 7?. if and only if d < d. This proves 
the statement of the corollary. □ 

In order to infer Theorem 3 from Corollary 6 we still need to remove the assumption that 
equalities be preserved along the evolution. To do that, we discuss the two cases in which this 
assumption could be restrictive, that is, when the solution starts at a discontinuity because 
either xb — %o = 1 or xq — xa = 1. 

i) We assume that xo(0) = xa(0) + 1 and xb — xq > n ^ j X ■ (Indeed, when xb — xa < 
1 + n ^ } 1 we already know that there is a solution leading to cluster coalescence.) Then, 
we note that io(0) = ub(xb — xq) — Ua while x a (0) € [0,1] for every a <G A. As 
xb — xq > "rc^j" 1 implies xo (0) > 1 , the solution may not stay on the discontinuity, and 
in particular is such that xo(t) — x a (t) > 1 for every t > 0. 

ii) We assume that xq(0) = xb{0) — 1- Then, we note that io(0) = tib — ha{xq — xa) while 
ih(0) G [—1,0] for every b G B. As > tia implies xq(0) > 0, the solution may not stay 
on the discontinuity, and in particular is such that xs{t) — xo(t) < 1 for every t > 0. 

In view of this discussion, we infer that it is not restrictive for the robustness analysis to 
assume that equalities be preserved, and we conclude that the state of the perturbing agent 
xo can be chosen in such a way to make the two clusters merge if and only if xb — xa < d. 
This proves the statement of Theorem 3 when tia < Ub- The case of ua = n B is much 
simpler and is left to the reader. 
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